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ABSTRACT 

We present detailed computations of the vertical structure of an accretion disc 
illuminated by hard X-ray radiation with the code TITAN-NOAR suitable for Compton 
thick media. The energy generated via accretion is dissipated partially in the cold disc 
as well as in the X-ray source. We study the differences between the case where the X- 
ray source is in the form of a lamp post above the accretion disc and the case of a heavy 
corona. We consider radiative heating via Comptonization together with heating via 
photo-absorption on numerous heavy elements as carbon, oxygen, silicon, iron. The 
transfer in lines is precisely calculated. A better description of the heating/cooling 
through the inclusion of line transfer, a correct description of the temperature in the 
deeper layers, a correct description of the entire disc vertical structure, as well as the 
study of the possible coronal pressure effect, constitute an improvement in comparison 
to previous works. 

We show that exact calculations of hydrostatic equilibrium and determination of 
the disc thickness has a crucial impact on the optical depth of the hot illuminated zone. 
We study the lamp post model for a low (to = 0.03) and high (m = 0.3) accretion rate. 
In both cases we assume a moderate illumination where the viscous flux equals the 
X-ray radiation flux. A highly ionized skin is created in the lamp post model, with the 
outgoing spectrum containing many emission lines and ionization edges in emission 
or absorption in the soft X-ray domain, as well as an iron line at ~ 7 keV consisting 
of a blend of low ionization line from the deepest layers and hydrogen and helium 
like resonance line from the upper layers, and almost no absorption edge, contrary to 
the case of a slab of constant density. A full heavy corona completely suppresses the 
highly ionized zone on the top of the accretion disc and in such case the spectrum is 
featureless. 

Key words: galaxies: active - accretion, accretion discs - black hole physics - galax- 
ies:Seyfert, quasars - X-rays. 



1 INTRODUCTION 

Broad band spectra of active galactic nuclei (AGN) demon- 
strate the presence of both cold and hot plasma in the vicin- 
ity of a black hole. Furthermore, detailed analysis of X-ray 
spectra reveal that those two phases interact radiatively dur- 
ing their accretion on a central object (for recent review, see 
Collin et al. 2000, Reynolds 2000; see also Poutanen 1999 
for the discussion of this problem in case of galactic black 
holes). However, the accretion pattern is still a subject of 
dispute and several detailed scenarios were proposed. 

In most cases, a key element of the model is an accre- 
tion disc which extends, or does not, down to the marginally 
stable orbit. The disc surface is expected to be strongly irra- 



diated by the X-ray emission containing a significant fraction 
of the observed bolometric luminosity of AGN. 

At a given disc radius, two distinct situations are possi- 
ble. The hot plasma may provide only a source of irradiating 
photons, being located well above the disc or at different ra- 
dius (the 'lamp post' model), or it may be located at the con- 
sidered radius, in hydrostatic equilibrium with the cold disc, 
additionally exerting a dynamical pressure on the cold disc 
surface (the 'heavy corona' model) and modifying strongly 
its ionization state (Rozanska et al. 1999). 

In both families of models the optically thick disc is 
illuminated by a hard X-ray flux. Those two cases, however, 
differ significantly with respect to the ionization state of the 
disc surface layers, as indicated by the study of Rozanska et 
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al. (1999). The density at the surface of a disc, irradiated 
through some non-locally produced X-ray flux, is zero while 
the density of the disc with a corona is quite high due to the 
coronal pressure. This may lead to significant differences in 
the disc radiation spectra, both in the continuum and the 
properties of the iron K a line. 

Careful approach to the modeling of the broad band disc 
spectra is essential in order to test both models, particularly 
in view of the first chance for true X-ray spectroscopy, com- 
ing with XMM and Chandra observations. Therefore, in this 
paper, we address the problem of the vertical structure and 
the spectrum of the radiation emitted by the illuminated 
disc in two extreme cases: with and without the dynamical 
effect of the corona. 

For that purpose we consider the radiation transfer 
within the disc much more carefully than it was done by 
Rozanska et al. (1999), applying to the surface layers the 
code titan developed by Dumont, Abrassart & Collin (2000 
hereafter DAC00). This code also provides us with the con- 
tinuum and line spectrum of the disc and is interfaced with 
a Monte Carlo code NOAR to take into account Compton 
diffusion properly, and particularly to compute the iron K a 
line profile and the high energy continuum. 

This approach itself is an improvement in comparison 
to the earlier results on the radiation spectra of X-ray ir- 
radiated accretion discs, although it is still not free from 
important simplifications. 

We first stress that all computations of disc structure 
and disc spectrum are (and will necessarily be) based on 
simplifying assumptions because the problem is too compli- 
cated to be solved exactly. And as it is not possible to make 
a 'uniform' improvement - it is worth to improve one aspect 
of the model at the expense of another one. But different 
approximations leading to different predictions, it is neces- 
sary to perform a detailed comparison of the results, and 
this is not possible unless exactly the same model is used in 
the computation. For instance Ballantyne et al. (2001) and 
Nayakshin et al. (2000) have computed the spectrum of an 
irradiated disc in hydrostatic equilibrium using codes simi- 
lar to ours. In the most simple case (an irradiated slab with 
a constant density), with exactly the same parameters, the 
reflected spectrum predicted by the three codes shows im- 
portant differences in the detailed spectral features as well 
as in the overall shape of the continuum (cf. Pequignot et 
al. 2001). Thus in order to judge the reliability of the results 
and to compare them to the observations it is absolutely 
necessary that all authors be very clear on the approxima- 
tions made in their computations, so that everyone could 
be aware of them. This is what we try to do in the present 
paper. 

As improvements with respect to other methods, we 
have chosen here mainly to favor the radiative transfer treat- 
ment (other codes use the escape probability approximation, 
especially for the lines), and to include a correct (but obvi- 
ously model-dependent) description of the vertical structure 
of the viscous disc underlying the irradiated layers (other 
codes use rough estimates from vertically averaged models) . 
These improvments are made at the expense of the atomic 
data. However though our computations include less atomic 
transitions than XSTAR (used by Nayakshin et al. 2000), 
they include more ions and transitions than Ballantyne et 
al. 2001. 



In our approach the viscosity in the disc is proportional 
to the total (gas + radiation) pressure. Close to the sur- 
face we compute the radiative transfer of continuum and 
lines using two-stream approximation with non-LTE treat- 
ment and with heavy elements abundances as described in 
DAC00. Deeper inside the disc the diffusive radiative trans- 
port is assumed. Due to better heating/cooling description 
we are able to calculate regions with temperatures below 
~ 10 5 , well matching with the disc interior, i.e. full profiles 
of the disc from the surface to the equatorial plane. 

We compute the emergent spectrum from the UV to the 
hard X-ray range, including features like the iron K a line. 

This paper is restricted to a moderate illumination (the 
lamp post and the heavy corona models), with low and high 
accretion rate, while the following paper will be devoted to 
the high illumination case (the flares). 

In section 2 we present the assumptions and numerical 
method of computation. The results are presented in section 
3. The final sections are devoted to discussion and conclu- 
sions. 



2 METHOD 
2.1 Generalities 

We consider the vertical structure of an illuminated accre- 
tion disc and the emitted spectra in two cases. In the first 
case hard X-ray emission is produced in an active region lo- 
cated above an accretion disc (the lamp post model), and in 
the second case X-rays come from a hot corona (the heavy 
corona model) . The spectral distribution of X-rays is a free 
parameter of our model and in both cases it is assumed to 
be the same. The only difference is that corona is heavy 
and exerts a considerable pressure on the disc surface thus 
modifying its properties. 

In the case of lamp post model the surface boundary 
conditions are the same as for a standard accretion disc, 
i.e. the density (and gas pressure) approaches zero at the 
surface. 

In the case of corona above an accretion disc the elec- 
tron temperature of the hot plasma is chosen arbitrarily, in 
agreement with observations. We do not consider any par- 
ticular coronal model in this paper. The optical depth and 
the gas pressure of the corona are determined according to 
the spectral distribution and the electron temperature from 
the relation that the flux emitted by the corona Fx is pro- 
duced via Compton cooling of hot electrons by the soft flux 
Fsoft outcoming from the accretion disc. The non-zero value 
of the coronal gas pressure implies surface condition on the 
disc atmosphere. 

The disc structure in both cases is calculated in the 
same way, and the presence or absence of the heavy corona 
influences the disc only through the surface boundary con- 
ditions. These computations require an assumption about 
the viscous energy generation within a disc. We adopt a 
standard a viscosity model where the viscous stress is pro- 
portional to the total pressure. 

The results depend on the mass of the central black 
hole, M, the distance from the center, r, the total flux gener- 
ated by accretion, F gen (or total accretion rate m), the frac- 
tion of energy generated in the X-ray source, / = Fx/F gen , 
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and its spectral shape, and finally on the viscosity parame- 
ter a. Since we do not specify any coronal model, we assume 
in both cases that all the angular momentum is transported 
outwards by the disc, but only a fraction of the gravitational 
energy, 1 — /, is dissipated in the disc. Therefore, the disc 
accretion rate is smaller than the total accretion rate which 
is the parameter of our model, because the flux dissipated 
in the disc, Fdisc, is smaller than the total generated flux; 
Fdisc = F gm — Fx ■ Note that in the case of lamp post model 
/ is not equal to Lx / Ltoi ■ 

The computations of the disc structure are done 
through iterations between the code solving the hydrostatic 
equilibrium equation and the code solving the radiation 
transfer, as described below. 



2.2 The computation scheme 

2.2.1 The beginning step 

The first step of the computations is completed as described 
by Rozanska et al. (1999 Section 2). The radiative transfer 
within the disc is treated in the diffusion approximation and 
the convective transport is also included. The viscous energy 
dissipation is given by the local a viscosity description. The 
hydrostatic equilibrium completes the set of equations which 
fully determines the vertical temperature and the density 
profiles. 

Thus, in the first step, we integrate the set of equations: 
„ _ 16<tT 3 dT _ ^ _ 
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where: F, F rad , F conv , T, p, z, P, P ga s, Prad, ftjf, a and r arc 
respectively: the flux emitted by the disc, the flux carried by 
radiation, the flux carried by convection, the temperature, 
the density, the vertical coordinate, the total pressure, the 
gas pressure, the radiation pressure, the Keplerian angular 
velocity, the Stefan-Boltzman constant and the Rosseland 
mean optical depth. 

The opacity k (the Rosseland mean) as a function of 
density and temperature is taken from Alexander, Johnson 
& Rypma (1994) for log T < 3.8, from Seaton et al. (1994) 
for log T > 4.0, and it is interpolated between these two 
tables for intermediate values of the temperature. 

For convection we adopt a simple description based on 
the mixing length theory used in stellar interiors. 

For both models we assume that half of the flux Fx is 
directed toward the disc and this fraction is absorbed (the 
albedo equals zero in the first computational step). There- 
fore, in both cases we have the same boundary conditions 
on the flux and temperature at the surface. Using the Ed- 
dington approximation we require: 



F(H d ) = F soft ee 0.5F X + (1 - f)F gen , 
oT 4 (H d ) = 2°" T e 4 // = 2 F sofU 



(6) 
(7) 



where Hd is half of the disc thickness actually determined 
by the boundary condition on the equatorial plane: 



F(z = 0) = 0. 



(8) 



The third boundary condition is different for both mod- 
els. In the lamp post model, the density at the disc surface 



p{H d ) = 0. 



(9) 



In the heavy corona model this condition results from 
the value of the pressure at the basis of the corona and the 
requirement of pressure equilibrium between the disc and 
the corona: 



p(Hd) = ft, 



prn H 
k B T{H d y 



(10) 



Assuming the vertical component of the gravity to be ver- 
tically constant, the hydrostatic equilibrium of the corona 
leads to the gas pressure of corona: 

-.2 T c 



cor — ^ K H-cori 



(11) 



where the pressure scale-height of the isothermal corona is 
of the order of: 

"—(St)""' 1121 

where the opacity for electron scattering, ks - the 

Boltzman constant and mn - the mass of hydrogen. We 
assume the value of the mean molecular weight p = 0.5 
appropriate for cosmic chemical composition. Note that the 
pressure scale height H cor would be lower than in the case of 
variable gravity by a factor of \[2. Those two previous equa- 
tions are equivalent, but they show clearly the dependence 



of ft. 



and H cor on the free parameter T cor . The coronal 



plasma is taken to be a single temperature medium with 
T cor equal to the electron temperature chosen arbitrarily at 
1 x 10 9 K (as suggested by observations of high energy cut off 
of X-ray spectrum). The optical thickness of the corona r c 
is not a free parameter of our model, but it is determined by 
the energy equilibrium. Since the corona is cooled via seeds 
photons coming from the accretion disc, it should have a 
specific thickness to produce an X-ray power-law spectrum 
with the energy index chosen in our computations. We con- 
sider the case of an energy index cue = 0.9, and we compute 
the coronal optical depth using a simple Comptonization 
code (Czerny & Zbyszewska 1991) based on semi-analytical 
formulae and appropriate for optically thin media. 

We take into account only gas pressure due to the 
corona in this first step. The radiation pressure due to il- 
lumination is included in the second step of iteration after 
solving radiative transfer. 

We solve this two point boundary problem by a shooting 
method. The integration of equations presented in this sec- 
tion is performed by the second order Runge-Kutta scheme 
with adaptive stepsize, as in Pojmanski (1986). 



2.2.2 Radiative transfer in the surface layers 

Knowing the vertical density profile we can solve the ra- 
diative transfer within the surface disc layers much more 
accurately using the code titan/noar of DAC00. titan is 
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based on Eddington two-stream approximation of the radia- 
tive transfer and works in plane-parallel geometry. Radiative 
transfer is computed both in lines and in continuum (i.e. 
lines are not treated with escape probability approach) so 
the code can be used for very inhomogeneous thick media, 
as in the present case. 

Thermal equilibrium and ionization and statistical equi- 
librium of ions are computed in complete non-LTE for 10 
elements and all corresponding ions. We consider following 
elemental abundances: H: 1, He: 0.085, C: 3.3 x 10~ 4 , N: 
9.1 x 10~ 5 , O: 6.6 x 10~ 4 , Ne: 8.3 x 10~ 5 , Mg: 2.6 x 10~ 5 , 
Si: 3.3 x 10" 5 , S: 1.6 x 10~ 5 , Fe: 3.2 x 10" 5 , from Allen 
(1973). Comptonization is taken into account both in the 
thermal equilibrium and in the computation of the spec- 
trum, through an iteration with the Monte Carlo code of 
Compton scattering, NOAR. 

The radiative transfer equations for each frequency v 
can be written in the form: 



1 dl+ 

VS dz 



-{k v + a es )Iv + 



1 dl v , \ | — I °es r -| 



(13) 



(14) 



where z is the distance measured from the disc surface, k u is 
the absorption coefficient in cm -1 , a es is the electron scat- 
tering coefficient, and r\ v is the emissivity. 

The boundary conditions are imposed at the illumi- 
nated side of the disc surface and asymptotically in the disc 
interior. 

The surface boundary condition is given by: 



/, + (0) 



2n " 



(15) 



where F„ is the frequency dependent illuminated X-ray 
flux. 

It is not necessary to continue with integrations down 
to the equatorial plane of the disc since non-LTE effects 
disappear completely when the optical depth is ~ 10 , and 
where the diffusion approximation is satisfactory anyway. 

Therefore we divide the disc arbitrarily into a surface 
layer of thickness d and the interior, like in the case of a 
stellar atmosphere problem. The disc interior is calculated 
in the diffusion approximation of radiative transfer, and the 
zone down to z — d is computed by titan. However, we have 
to impose a second boundary condition at z = d which will 
insure that the energy flux dissipated in the disc F^ tsc {d) 
will leave the interior: 

I-(d) = I+{d) + ^-Ff sc (d). (16) 

Since the disc computations provide us only with the 
frequency-integrated value of the flux Fduc dissipated in the 
disc and the value of the temperature at z = d, we determine 
the spectral shape of the flux using the non-grey diffusion 
approximation but neglecting the frequency dependence of 
the opacity coefficient 



F?*°(d) = f t^p-=F disa {d) 



■kB v (T) xexp(cc) 
4<r 4 (exp(z) - 1) 



(17) 



where x = hv/ksT, B V (T) is the Planck function. This con- 
dition is equivalent to the condition used by NKK00 (Eq.14) 
in their illuminated disc model. 



The code returns the temperature profile and the spec- 
tral shape of the emitted radiation. The results are not sen- 
sitive to the choice of d as long as T es (d) is neither too small 
(smaller than 3) nor too large (larger than 30) since in the 
first case the radiative transfer of the incident X-ray flux is 
not described accurately, and in the second case the accu- 
mulated error of the computations made practically in LTE 
zone is too large. We arbitrarily choose r es (d) ~ 6. In ear- 
lier works it was estimated that X-rays penetrate the disc 
atmosphere up to an optical depth approximately equal 3 
(Sincell & Krolik 1997, NKK00) 

2.2.3 Density profile from hydrostatic equilibrium 

Having the new temperature profile of the surface layers 
from the code titan calculated for the density profile in 
the first step, we can calculate a new density profile from 
hydrostatic equilibrium. For that purpose we solve the e ntire 
disc structure basically as described in Section 2.2.1, but 
with few modifications. At a distance smaller than d from 
the disc s urface we assume the temperature profile given 
by Section 2.2.2] and we do not solve the diffusive radiative 
transfer (Eqs. 1-3), while deeper within the disc we solve 
the full disc vertical structure equations including diffusive 
radiative transfer and heat generation. At this stage we also 
include the radiation pressure due to the incident radiation. 
Having the radiation field given by TITAN we determine this 
pressure component by: 

dP r ad 



dz 



1 

cV3 



k v F v dv, 



(18) 



where the radiation pressure on the surface is equal to Fx / c. 

In the heavy corona model the value of the density at 
the disc surface changes in subsequent iterations due to the 
change of the surface temperature, since it is determined 
from the adopted value of the pressure at the basis of the 
corona. This differs from the case of the lamp post model 
where the surface density is zero through all iterations. 

One full iteration is done after computing the tempera- 
ture and density profiles. The new density profile is used to 
repe at rad iative transfer computations as described in Sec- 
tion 2.2.2. We again divide the disc in a surface layer with 
optical thickness r es (d) ~ 6, and the rest of the disc, where 
the diffusion approximation is adopted. 

We repeat iterations between radiative transfer and 
density profile until the temperature profile does not change, 
and in section below we present the results. 

2.2.4 Iron li- ne 

When the irradiated disc structure is determined we cal- 
culate in more detail the shape of the X-ray spectrum in 
order to determine more accurately the properties of the 
iron K a line. These computations are done by the Monte 
Carlo code NOAR described in detail by DAC00. This code 
describes the radiative transfer of line components coming 
from various iron ions, as determined by titan, and includes 
the line broadening due to Comptonization. It is used also 
after each titan run to provide the net (heating - cooling) 
Comptonization effect. We do not include here any kine- 
matic broadening related do the systematic motion of emit- 
ting material thus providing an 'intrinsic spectrum' which 
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Figure 1. Temperature profile versus Thomson optical depth for 
the lamp post model. The initial profile is denoted by the dotted 
line, next iterations are presented as long dashed lines. The final 
profile is shown by the solid line. The vertical straight line marks 
the transition between the diffusion approximation and the exact 
solution. 
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Figure 2. Temperature profile versus Thomson optical depth for 
the lamp post model. The initial profile is denoted by the dot- 
ted line and the final profile is shown by solid line. The vertical 
straight line marks the transition between the diffusion approxi- 
mation and the exact solution. 



later may be folded with arbitrary flow pattern like Keple- 
rian disc motion, inflow or outflow. 



3 NUMERICAL RESULTS 

All results are presented for a single representative radius 
r — WRschw (where Rschw = 2GM/c 2 ), for a mass of black 
hole M = 10 s Mq, two accretion rates (in units of Eddington 
accretion rate with efficiency 1/12): rh = 0.03 and rh = 0.3, 
a viscosity parameter a = 0.1, and for / = 0.5. For these 
parameters, the total X-ray flux as well as locally dissipated 
flux are equal to 6.96 x 10 13 erg s 1 cm" 2 , with an incident 
X-ray radiation flux equal to half of this value. 

The distribution of incident X-ray spectrum is a power- 
law with energy index oe = 0.9, extending from 2.8 eV up 
to 100 keV. 

3.1 Vertical structure of the disc without corona 

We perform the coupled computation of the hydrostatic 
equilibrium and the radiative transfer for an irradiated ac- 
cretion disc as described in Sect. ^. 

We do not include conduction but the iterative method 
allows to pass effectively from the hot upper layers to the 
cold inner disc without serious instability problems. Con- 
trary to the computations done by NKK00 we do not start 
calculation of radiative transfer from an already unstable 
temperature profile determined with the photoionization 
code. The successive iterations of the temperature profile 
for rh = 0.03 are shown in Fig. |l| starting from the solution 
given by diffusion approximation (dotted line). 

The temperature distribution displays an expected pat- 
tern; the hot upper layer reaches approximately the inverse 
Compton temperature (determined by the spectrum of the 
incident radiation and the radiation spectrum emitted by 
the disk and its atmosphere) , then follows a rapid tempera- 
ture drop due to an increase of the density and, consequently, 



an increase of the cooling efficiency of the gas. The deeper 
layers of the disc are not affected by the irradiation, so the 
better description of the radiative transfer than in the pa- 
per of Rozariska et al. (1999) does not alter their results for 
the disc interior. After an initial decrease close to the sur- 
face, the temperature inside the disc is rising towards the 
equatorial plane, reaching T ~ 2.94 x 10 5 K. The solution 
converges satisfactorily after seven iterations. The number 
of iterations cannot be easily increased as the computations 
are time-consuming, taking about 3 days per iteration. 

The temperature profile for the case of high accretion 
rate is presented in Fig. ^| The heated zone is thicker in this 
case than in the lower accretion rate solution. 

In the final stage of iteration with radiative transfer we 
obtain an optical depth of the Compton heated skin t^°* ~ 
0.03 (defined by a decrease of the temperature by a factor 2 
with respect to the surface temperature, in reference to the 
simple approach of Begelman, McKee and Shields 1983, the 
discussion of McKee and Begelman 1990 and Appendix A) 
for rh = 0.03 and r e h s ot ~ 0.048 for rh = 0.3 (see Figs. | and 
[|. For m = 0.03 the height of the bare disc is H d ~ 3.7x 10 12 
cm, and the height of the hot skin is H^ot ~ 5.1 X 10 12 cm. 
The transition layer is very thin, with Htran — 6 x 10 11 
cm. The same geometrical sizes in case of m = 0.3 are as 
follows: H d ~ 2.3 x 10 13 , H hot ~ 1.44 x 10 12 , and H tran = 
1.1 x 10 11 cm. These values are much smaller than the disc 
radius (3 x 10 14 cm) which means that both the disc and 
the hot layer are geometrically thin. 

Equation |l^ gives analytically the scale-height of the 
Compton heated zone, when T cor is replaced by the inverse 
Compton temperature Tic- For the accretion rate rh = 0.03 
and Tic — 2.1 x 10 7 K, as calculated from the spectral 
distribution of hard X-rays and the disc radiation, we can 
analytically derive the value of H^ot = 2.6 x 10 12 cm. We 
see that the actual thickness of the Compton skin is a factor 
of 2 larger when exact vertical calculations are done. The 
geometrical extension of the hot zone formally depends on 
the adopted value of the initial small (but non-zero) density 
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Figure 3. Density distribution versus Thomson optical depth for 
the lamp post model. The initial profile is denoted by dotted line, 
next iterations are presented as long dashed lines. Final profile is 
shown by the solid line. 
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Figure 4. The pressure distribution versus Thomson optical 
depth in the surface layers for lamp post model after the last 
iteration for m = 0.03. 



at the surface but in practice this dependence is weak since 
the gravity in the disc increases outwards. 

In Fig. ^ we show the density profile for rh = 0.03, again 
starting from the diffusion approximation. The density of 
the illuminated layer is lower than for the bare disc up to 
a Thomson optical depth ~ 0.1. Deeper in the disc, the 
density profile causes a density inversion due to convection, 
as shown in Rozanska et al. (1999). 

The final gas pressure and total pressure profiles for 
m = 0.03 are shown for surface zones in Fig. ^. Radiation 
pressure always dominates the gas pressure by at least two 
orders of magnitude. 

The pressure profile resulting from diffusion approxi- 
mation goes monotonically through the disc. Nevertheless 
in outer layers we see a 'wiggle' in the region of the temper- 
ature drop but the hydrostatic equilibrium and the radiative 
transfer are nevertheless computed with the required accu- 
racy at each step. We argue here, that such 'wiggles' are due 
to thermally unstable regions, which are expected to appear 
in an atmosphere in hydrostatic equilibrium and illuminated 
by hard X-rays (RC96). 



Since we start our computation from the thermally sta- 
ble temperature profile of the non illuminated disc atmo- 
sphere, and since we perform an iteration between this pro- 
file and a single value density profile, we do not achieve a 
three value temperature profile (i.e. thermally unstable re- 
gions) after the radiative transfer computation. Note that 
iterating between temperature and density profiles is a usual 
way of treating any radiative transfer computations in hy- 
drostatic equilibrium (see Hubeny 1990, Madej & Rozanska 
2000, NKK00). In thermally unstable regions, pressure is a 
monotonic function of optical depth (RC96), and therefore if 
we keep temperature as a single value function of r es in our 
radiative transfer computations, the pressure inversion will 
appear in the thermally unstable region. We note, however, 
that the wiggle occupies only a small part of the irradiated 
layer and in particular it does not affect the overall shape of 
the emitted spectrum but it may h ave s ome impact on the 
detailed spectral features (cf. Sect. 4.1.3). 



As a conclusion we stress that in any radiative trans- 
fer computation of disc atmospheres illuminated by hard 
X-rays, if the initial temperature profile does not displays 
instabilities, the final pressure profile will display inversions. 
It may be the way to check that the cooling and heating are 
correctly computed, as the instabilities are due to discon- 
tinuities in the cooling rate through the illuminated atmo- 
sphere. Including the conduction flux into the scheme would 
solve the problem but it is technically difficult and no ex- 
act radiative transfer computations were performed with the 
conduction effect included up to now. 



3.2 Outgoing spectrum in the lamp post model 

3.2.1 Overall spectral distribution 

The reflected spectrum for rh = 0.03, obtained from the 
final iteration, is presented as a thick solid line in Fig. ^, 
which gives also the incident illuminating spectrum (power 
law line). The dashed line shows the observed spectrum, 
equal to the sum of the incident and the reflected spectra. 
Plotted spectra are usually degraded to a spectral resolution 
of R = 30, and no kinematic corrections are applied. 

The reflected spectrum consists of a large opti- 
cal/UV/soft X-ray component due to the disc emission as 
well as to thermalization of a significant fraction of the in- 
cident radiation. The X-ray part of the spectrum shows the 
characteristic spectral shape of reflection from partially ion- 
ized medium, with much more radiation in the soft X-ray 
band than in the case of a neutral medium (Lightman & 
White 1988), and in particular a lot of emission lines and 
ionization edges. 

The reflected spectrum is very different from that of 
a constant density slab. It is actually not easy to compare 
precisely the two cases, as for hydrostatic equilibrium the 
ionization parameter has no meaning and it is also difficult 
to define an average density. To perform the comparison, 
we have used a slab of constant density equal to 10 12 cm -3 
with exactly the same irradiation as in the hydrostatic case: 
a power law spectrum with the same flux incident on the top 
of the layer, and a thermal spectrum incident on the bottom 
of the layer, with F^ sc given by Eq. [l7| This is comparable 
to the method used for instance by Ross et al (1993) and 
Zycki et al (1994) to compute the disc spectrum, except 
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Figure 5. The reflected spectra of the irradiated disc for the lamp post model for hydrostatic equilibrium (thick solid line), and for 
constant density model (thin solid line). Final spectra are presented by dotted and dashed lines respectively. Power law line marks 
spectral distribution of incident radiation. Spectral resolution: 30. 



that they had a larger density, representative of the mid- 
plane value, while we have a density more representative 
of the irradiated layers. The resulting reflected spectrum 
is displayed as a thin solid line on Fig. [|, in order to be 
compared with the lamp post reflected spectrum, and the 
corresponding observed spectrum is shown as a thin dotted 
line. 

We see that the constant density case shows a larger 
thermal UV bump, very close to the underlying black body 
continuum, with only small traces of ionization edges in ab- 
sorption, while the hydrostatic model exhibits a Lyman edge 
in absorption and a strong absorption above 20 eV, caused 
by the presence at 0.1 r es ^ 1 of several weakly ionized 
species: O n and O III, S II and S in, C in, Si II and Si ill 
ions. As an illustration Figs. ^| a and b show the fractional 
abundances of the different iron ions as functions of the op- 
tical depth, respectively for the hydrostatic equilibrium and 
for the constant density case: in the latter case we see that 
only relatively highly ionized species are present in the deep 
layers, contrary to the hydrostatic case. In both cases highly 
ionized species are present at r es ^0.1. 

The low ionization state in the hydrostatic case is due 
to the high density of the deep layers. Since LTE is almost 



reached, the ionization stage depends only on the density 
and the temperature profile. It is therefore most important 
to solve the whole vertical structure of the disc to determine 
the abundances of these weakly ionized species and to get the 
correct spectral distribution in the far UV range. 

In the soft X-ray range the shape of the reflected spectra 
differs also in the two cases, as it is steeper for the constant 
density. This is due to the lack of thermal instability and 
consequently the smooth distribution of the temperature in 
the constant density case. On the contrary, the shape of the 
reflected spectrum in the hard X-ray range is almost the 
same in the two cases, as the Thomson thickness of the hot 
layer is relatively small (< 0.1) so Compton reflection takes 
place in a cold medium. Also because of the small optical 
thickness of the hot skin, the outgoing spectrum is not mod- 
ified by Comptonization in the UV and soft X-ray ranges. 
Finally, one should notice that in the soft and hard X-ray 
range the overall shape of the observed continuum is quite 
similar in the two cases, and it is only the (very different) 
detailed spectral features which could help to distinguish be- 
tween them. In particular, the spectrum in the hydrostatic 
case exhibits strong ionization edges in emission in the soft 
X-ray range, but no ionization edge of Iron in absorption 
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Figure 6. The fraction of iron ions in case of lamp post model 
versus Thomson optical depth in constant density case (upper 
panel) and hydrostatic equilibrium (lower panel). Accretion rate: 
rh = 0.03. 



at ~ 10 keV, contrary to the constant density case (see 
Figs. |E] and ^); one can see only a weak Fexxvi ionization 
edge in emission. Computations performed under assump- 
tion of constant pressure give results intermediate between 
the constant density and hydrostatic equilibrium solutions 
(Dumont et al. 2001). 

The effect of the adopted value of the accretion rate 
can be seen comparing Figs, [s] and |?]. The spectra are qual- 
itatively similar, but they also show systematic differences 
matching the difference in the temperature profile. 

The overall reflectivity in the X-ray band is higher for 
higher accretion rate, as it corresponds to somewhat higher 
ionization state of the gas. This is well seen in the plot of 
fractional abundances of iron (see Fig. |s|). FeXXV is emitted 
predominantly from the optical depth of ~ 0.09 for rh = 0.03 
but from the optical depth of ~ 0.18 for rh — 0.3 while FeV 
comes from r es > 0.3 and r es > 0.7, correspondingly. There 
is systematically even more reflection in the soft X-ray band, 
making this case more similar to the constant density case 
than it was for lower accretion rate. Also for rh = 0.3 there 
are copious soft X-ray lines from CNO elements and other 
species, although they are generally less intense than for 
rh — 0.03 if measured with respect to the reflected com- 
ponent but have similar EW measured with respect to the 
total continuum. Higher average albedo leads to the reduced 
dominance of the optical/UV bump. 




Figure 7. The reflected spectrum (continuous line) and final 
spectrum (dashed line) for the lamp post model in hydrostatic 
equilibrium for high accretion rate. Spectral resolution: 30. 
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Figure 8. The fraction of iron 


ions in case 


of lamp post model 



versus Thomson optical depth in hydrostatic equilibrium. Accre- 
tion rate: rh = 0.3. 



3.2.2 K a line properties 

The differences of ionization state in the two cases induce 
quite different Iron spectrum, as already stressed by NKK00 
(who have performed a comparison using a constant irradia- 
tion flux and not a constant Fx /F V i SC ratio like here). Figs. ^ 
and [l^ display enlarged versions of the reflected spectra be- 
tween 3 and 15 keV with relatively high spectral resolution. 
We see that in the hydrostatic case the spectrum is domi- 
nated by two components: the 6.4 keV component is due to 
several Iron Ka lines of Fexm and less ionized species (the 
most intense line being due to Fev), and the 6.7 keV compo- 
nent is the Fexxv recombination line blend with Fexxiv. For 
rh = 0.03 the first one slightly dominates while for rh = 0.3 
the second one is stronger. The recombination line of Fexxvi 
is also present at 6.9 keV but is much weaker, as well as the 
K/3 line at 7.1 keV. An emission edge due to Fexxv is seen at 
8.8 keV. It is difficult to compute the equivalent width EW 
of each component separately, as they are blended together 
through Compton broadening. The whole line complex has 
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Figure 9. The close-up of the iron line region of the reflected 
spectrum (in arbitrary units) for lamp post geometry and con- 
stant density case (dotted line) and hydrostatic equilibrium (solid 
line). Spectral resolution R = 100. Accretion rate: m = 0.03. 
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Figure 11. The temperature distribution for the heavy corona 
model. Initial profile is represented by dotted line, and final iter- 
ation by solid line. 
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of the FeXXVI line is overestimated, as it integrates a blue 
wing partly due to inverse Compton broadening of the other 
lines, and as a consequence the other lines are slightly un- 
derestimated. 

In the constant density case the spectrum is dominated 
by the Fexxv line and by a corresponding intense ioniza- 
tion edge in emission. There is also a smooth but strong 
absorption above 8 keV due to the superposition of several 
ionization states. 
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Figure 10. The close-up of the iron line region of the reflected 
spectrum (in arbitrary units) for lamp post geometry and hydro- 
static equilibrium (solid line). Spectral resolution R = 100. 

an EW of 550 eV with respect to the reflected continuum, 
corresponding to an EW of 90 eV with respect to the ob- 
served continuum in the case of the lower accretion rate. For 
high accretion rate, the corresponding numbers are 405 eV 
and 90 eV. Dividing roughly the whole line profile into the 
three lines, we obtain the values of the EW of the line com- 
ponents given with respect to reflected and total observed 
continuum, correspondingly (see Table [j]) . 

These numbers are approximate: in particular the EW 



Table 1. Iron K a line properties in the considered models*. 
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3.3 The weight of the corona effect 

3. 3. 1 Vertical structure of the disc with coronal influence 

In the case of the heavy corona we assume the electron tem- 
perature of the hot plasma to be equal 1 x 10 9 K as expected 
from observations. This temperature and the assumed spec- 
tral distribution of X-rays imply an optical depth of corona 
Tcor = 0.355 due to the role of Comptonization (Haardt & 
Maraschi 1991). Therefore, the coronal pressure, calculated 
from Eq. 0, is equal to 8200 in cgs units. All other param- 
eters are the same as described in Sect. 3.1 

The dynamical pressure of the corona modifies the 
disc atmosphere, making the gas relatively cool and dense 
(Rozariska et al. 1999). 

The uppermost layers are heated up only by a fac- 
tor of ~ 1.2 (Fig. 0). The weight of the corona modifies 
the density of the surface layers and we show this effect in 
Fig. [l2| The initial high value of the density decreases to- 
wards the disc interior forming the outer density inversion 
(see Rozariska et al. 1999). Such inversion is reflected in the 
pressure variation versus Thomson optical depth (Fig. |l2[ ). 
Deeper gas density inversion is again connected with con- 
vection. 

In the outer layers the gas pressure is comparable to 
the radiation pressure, but towards the equatorial plane the 
radiation pressure dominates even by two orders of magni- 
tude. 
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Figure 12. The distribution of density, gas pressure and total 
pressure for the heavy corona model. Initial profiles are presented 
by dotted line, and final iteration by solid line. 
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Figure 14. The fraction of abundances of iron ions for the heavy 
corona model, versus Thomson optical depth. 
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Figure 13. The spectrum of the irradiated disc with corona. 
Reflected spectrum - solid line, Comptonized reflected component 
- dashed line. Accretion rate: m = 0.03. 



Figure 15. The close-up of the reflected spectrum (in arbitrary 
units) in the energy band of the iron line, for the heavy corona 
model. Spectral resolution R = 100. Accretion rate: m = 0.03. 



3.3.2 Overall spectrum 

The final spectrum as seen by the observer has to be calcu- 
lated taking into account the direct presence of the corona. 

In the case of a lamp post model we had to include, in 
agreement with equation pi half of the incident X-ray radia- 
tion which goes directly to the observer (see final spectrum 
as a bold dotted line in Fig. |E|). In the case of the corona 
model the radiation emerged from the disc is subsequently 
Comptonized by the hot coronal layer. To describe this ef- 
fect we apply a simple code of Czerny & Zbyszewska (1991) 
based on semi-analytical formulae appropriate for optically 
thin Comptonized medium. This code does not include the 
effect of the anisotropy of Compton scattering. However, 
since in a bare disc model we also simplify the description 
of the fraction of the primary radiation going directly to 
an observer (we use the power law instead of Comptonized 
emission shape) such an approximation provides a better 
comparison of the two cases. 

In Fig. [l^ we present the reflected local spectrum by 
continuous line, while the dashed line represents reflected 



spectrum modified by its passing through the corona with 
temperature 1 x 10 9 K and with r cor = 0.355. 

Such a corona, however, cannot cover the entire disc sur- 
face. According to our model, the sum of the energy emitted 
in X-rays and the energy emitted by disc cannot exceed the 
total energy generated due to accretion. Taking the assumed 
fraction of energy dissipated in the X-ray source / = 0.5 and 
using Eq. ^ (albedo equal to zero) to describe the total out- 
going soft flux, the Comptonized flux is given by: 

Fx = A c (0.5Fx + 0.5F gen ) = A C (1.5F X ), (19) 

which gives the approximate constraint on the amplification 
factor: 

Ac = 0.667. (20) 

Such Compton amplification in case of plane parallel 
coronal slab would give a very steep soft spectrum, in- 
consistent with the slope assumed in the computations. 
Our Comptonization code indicates that for the parameters 
Tcor = 0.355 and T = 1 X 10 9 K (appropriate to reproduce 
the spectral index a = 0.9) the amplification is 3.9. Also 
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from simulations done by Abrassart & Czerny (2000, their 
Fig. 3) one can see that in case of a disc black body spectral 
distribution with the temperature 50 eV the amplification 
factor never drops below 4 for a wide range of coronal tem- 
peratures. This is the main reason why a uniform non-patchy 
corona covering the entire disc is not a good model of Seyfert 
galaxies. 

Since we model the entire spectrum at a single radius, 
we can circumvent this problem by breaking down the as- 
sumption of full, continuous corona and consider the case 
of a corona covering the disc only partially. We introduce 
a new parameter, f ge0 m, which decreases the probability of 
scattering of the soft photons in the corona due to its patchy 
structure. This parameter has to be included in Eq. as: 

Fx = Acf a e Om (0.5F x + 0.5F gen ) = A c f gBom {1.5F x ) , (21) 
and implies the amplification factor to be: 
0.667 



Ac 



fgeom 



(22) 



Varying f ge om we can find the solution which preserves the 
required spectral slope, and the amount of coronal material 
which has to cover the cold disc to maintain energy balance 
between the hot and cold phases. 

Fig. [u^] shows that the spectrum reflected from the disc 
part covered by the corona and the observed spectrum from 
that radius is clearly a reflection from a basically neutral 
gas. No strong soft X-ray lines, characteristic of the lamp 
post geometry, are visible. Even at the disc surface the iron 
is only weakly ionized (Fig. |l4"| ). Only the 6.4 keV iron line 
is visible in the expanded picture of the 4-10 keV band, with 
relatively strong line and considerable edge (Fig. [IS). 

However, the final spectrum of the object is not just a 
spectrum at the radius covered by the corona, but also the 
contribution from the uncovered part has to be included, in 
agreement with the required covering fraction. In our specific 
case with parameters: spectral slope a = 0.9 and fraction 
of energy dissipated in the corona equal 1/2, fgeom is very 
small, ~ 0.13, and the additional contribution (~ 0.87) from 
disc radii uncovered by the corona will dominate the total 
spectrum. In this case it would be impossible to distinguish 
between the lamp post and the patchy corona models. 

However, in sources with a weaker Big Blue Bump (i.e. 
a small fraction of energy dissipated in the disc) f ae0 m is 
about 0.5 and the spectrum from a heavy corona model 
would approximately look like the spectrum in the lamp post 
geometry but with soft X-ray lines suppressed by a factor 
of order of 2. In sources with softer (steeper) X-ray spectra, 
slab amplification is significantly lower and the dumpiness of 
the corona is not needed. Therefore, in soft spectra sources, 
whenever the disc is covered by the corona, we would not 
expect any soft X-ray lines. 

3.4 Comparison with observations 

The are several interesting new observations of iron line 
done by satellites as BeppoSAX, XMM and CHANDRA. 
They usually report that iron line is rather narrow with- 
out relativistic broadening (Yaqoob, George, Nandra et al. 
2001). Reeves et al. (2001) showed observations of the ra- 
dio quiet quasar Markarian 205 where two components of 
iron line are present: a narrow neutral line at 6.4 keV and 



a broadened line centered at 6.7 keV. They concluded that 
those observations are consistent with a disc origin only if 
the disc ionization is high enough to produce He-like iron, 
and if the narrow neutral component is produced elsewhere. 
But those conclusions were obtained by fitting the data with 
constant density models. Our model of illuminated disc in 
hydrostatic equilibrium shows that both components of iron 
line are present. The ionization state varies across the atmo- 
sphere and therefore we get a recombination iron line from 
uppermost layers and a fluorescence line from deeper zones. 
Also the EW of fluorescence line in our model (EW = 40 eV 
for m = 0.03) is not in bad agreement with the observations 
(see Reeves et al. 2001 Table 1) where EW is 46 eV. 

Another interesting observation, which could be fitted 
by our model, was reported by Comastri, Stirpe, Vignali et 
al. (2001). They observed the bright Narrow- Line Seyfert 1 
galaxy Ark 564 and they reported a strong iron line from He- 
like iron line at ~ 6.76 keV. The EW of the line estimated 
from an ionized reflection model is ~ 87 eV. In our model 
with high accretion rate we obtain predominantly 6.7 keV 
component, with EW of ~ 90 eV. However, we must notice 
that this source has much softer X-ray spectrum (photon 
index F ~ 2.5, Comastri et al. 2001, Turner et al. 2001). 

Using our model, the emission lines from other elements 
such as oxygen, nitrogen and carbon can be determined. 
Branduardi-Raymont et al (2001) have observed the spec- 
trum of MCG -6 -30 -15 in the soft X-ray range with XMM. 
They claimed that the very intense and broad features in 
this object cannot be attributed to a warm absorber, but 
are emission lines produced by an irradiated accretion disc 
and relativistically broadened. This result was not confirmed 
by the analysis of the CHANDRA data, as reported by Lee 
et al. (2001). Our model shows that irradiated accretion 
discs with moderate illumination cannot produce very strong 
emission lines. The results presented in this paper (Figs. ^ 
and Q) show that the most prominent line (doublet Nejc 
Fexvn) has EW ~ 18 eV (for both m = 0.03 and rh = 0.3) 
in comparison to the emitted continuum. 



4 DISCUSSION 

4.1 Precision of the method, advantages and 
drawbacks 

4-1.1 Radiative transfer 

In the present paper we do not recall the influence of the 
transfer treatment, which was discussed in DAC 2000 and 
in Dumont & Collin 2001. It was for instance shown that it 
is most important to take into account correctly the transfer 
of the diffuse continuum (which is not performed in Ballan- 
tyne et al. 2001), and that the escape probability approxi- 
mation gives incorrect results for thick inhomogeneous slabs. 
Not only this approximation can lead to differences in the 
line intensities, but also to differences in the energy balance, 
which is dominated by the line heating/cooling in a large 
fraction of the illuminated layer. The importance of these 
processes will be discussed in more detail in a forthcoming 
paper devoted to this issue (Coupe et al. 2002). Their correct 
description is clearly an advantage of our model. 
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Figure 16. Influence of the atomic data on the reflection spec- 
trum for the lamp post model with rh = 0.03. Solid line: multi 
level Li-like, He-like, and H-like ions, dashed line: only multi level 
H-like ions. Significant difference is only seen below 200 eV. 



4-1.2 Atomic data 

The importance of heavy elements in illuminated disc 
atmospheres was clearly shown by Madej & Rozanska 
(2000). Thermal absorption is completely inefficient in pure 
hydrogen- helium atmospheres, and incident X-ray radiation 
is redistributed by Comptonization. Iron rich models con- 
sidered by Madej & Rozanska (2000) include only bound- 
free absorption on iron, but they show completely different 
spectra with smaller UV bump and a lot of ionization edges 
which do not exist when Iron is not taken into account. 

Obviously the atomic models and atomic data deter- 
mine directly the line intensities. However they have also 
an influence on the fractional ionic abundances, through the 
ionization from, and the recombination onto, the excited 
states. Indeed close to LTE, i.e. in the dense layers of the 
atmosphere, ionizations from excited levels almost balance 
recombinations onto these levels, while in the hot skin, only 
recombinations onto excited levels are important. Improved 
atomic data have been recently implemented in the code 
(cf. Coupe et al. 2002). In particular all Li-like and He-like 
ions are now treated in a complete interlocked way with 
5 to 9 levels plus the continuum, while in the old version 
of titan only H-like atoms were treated with several inter- 
locked excited levels. This is still less than in XSTAR, used 
by NKK00, but much more than Ballantyne et al. (2001). 
These new data have been used to compute the spectrum 
in the lamp post moderate accretion rate case, and they are 
shown to have an impact on the spectral features in the 
EUV, but basically not on the vertical disc structure and on 
the overall continuum spectral distribution (cf. Fig. |l^) . 



4-1.3 Thermal instability issue 

Another problem in which one has to choose an approxima- 
tion is in the approach of the thermal instability. 

The vertical structure of the surface irradiated layers of 
accretion discs is very complex (Rozanska & Czerny 1996, 
hereafter RC96). Illumination of gas in hydrostatic equilib- 
rium by hard X-rays implies the presence of a thermal insta- 
bility. More precisely, in certain pressure range, the gas can 



be in three states of thermal equilibria, with an intermediate 
one being unstable. The problem is best seen as the presence 
of S-shape feature in the temperature versus the ionization 
parameter E plot (Krolik, McKee & Tarter 1981). There is 
no unique solutions for the temperature and density pro- 
file. Physically, we expect a transition from the upper stable 
Compton-heated branch to the lower stable branch deter- 
mined by atomic processes. This transition must happen 
somewhere in the multi-solution zone and it may be sharp 
or a two-phase medium may be present. 

Unique strict solution of the vertical structure can only 
be found if electron conduction, in addition to the radia- 
tive transfer, is included (Rozanska 1999). However, con- 
duction complicates the numerical scheme considerably and 
at present it was never combined with computations of ra- 
diative transfer. In Rozanska (1999) radiative transfer was 
replaced by the energy balance equation, where heating and 
cooling functions were determined by means of photoion- 
ization code CLOUDY appropriate only for moderately thick 
media (DAC00). 

Other papers neglected the conduction term but con- 
tained more advanced radiative transfer computations. 
These papers can be divided into two families, depending 
on the way of iterating between hydrostatic equilibrium and 
radiative transfer, which implies their approach to the ther- 
mal instability problem. 

First family solves strictly both the radiative transfer 
and hydrostatic equilibrium, iterating density while solving 
radiative transfer. In this case the problem of appearing mul- 
tiple solutions is solved by an arbitrary choice of the transi- 
tion pressure and introducing a discontinuity. This discon- 
tinuity, chosen in the first iteration, is preserved in later 
computations. This approach was taken by Ko & Kallman 
1994, and subsequently by Sincell & Krolik (1997, 1999) and 
Nayakshin, Kazanas & Kallman (2000, hereafter NKK00). 
The choice of the position of this discontinuity does not in- 
fluence strongly the resulting spectra (Jimenez-Garate et al. 
2001). 

Second family keeps the density constant while iterat- 
ing the temperature at a given computation step. Such a 
numerical scheme always produces a unique solution across 
the transition region but at the expense of not providing 
a strict solution for both radiative transfer and hydrostatic 
equilibrium in illuminated atmosphere. This approach was 
used by Raymond (1993), and subsequently by Shimura, 
Mineshige & Takahara (1995), Madej & Rozanska (2000), 
Kawaguchi, Shimura & Mineshige (2001), and Ballantyne, 
Ross & Fabian (2001). 

Both approaches provide an approximate answer to the 
problem since the first one does not effectively iterate the 
optical depth of the hot zone adopted during the first iter- 
ation while the second one is never fully converged in the 
transition zone. 

In the present paper we have adopted the second ap- 
proach. This method gives a solution for the equilibrium 
temperature which corresponds to an ill defined value of the 
pressure. The extension of this layer is indicated in Fig. ^ 
which is an expanded version of Fig. [I], together with the 
extension of the multiple solution region. We see that the 
ill-defined region is relatively thin, so one would expect that 
it has no influence on the overall emitted spectral distribu- 
tion, but still it may lead to uncertainties in the detailed 
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Figure 17. The temperature profile for the lamp post model 
computed with rh = 0.03. Grey region: the multiple solution re- 
gion where the hot solution is chosen; hatched region: ill defined 
solution (see text). 

spectral features, anyway inherent to the presence of the 
thermal instability; in particular, it contains several impor- 
tant iron ions such as FeXVII. On the other hand the hot 
solution is suppressed at r es = 0.21, whilst the transition 
to the cold zone should probably occur at a slightly smaller 
optical thickness, when conduction is taken into account (cf. 
below). We have also checked that increasing the number of 
layers by a factor 2 in the transition zone does not change 
at all its structure. 



4-1-4 Expected effect of conduction 

Since our current approach does not include yet the conduc- 
tion term we test the possible effect of such a term using 
the local approach of Rozanska (1999) and computing the 
solution for rh = 0.03 with and without thermal conduction. 

Apart from the surface temperature given by the inverse 
Compton value, one more boundary condition was required 
to solve the second order differential equation of the energy 
balance with conduction. This condition was chosen as the 
requirement that at the bottom of the computed zone the 
solution matched that of the disc: when the same pressure as 
in the previous computations was reached, the temperature 
was also the same. 

The dashed line in Fig. ^ represents the temperature 
profile obtained without radiative transfer but from cooling 
and heating functions determined from CLOUDY (for bet- 
ter description see RC96). The solid line shows the same 
computations taking into account thermal conduction. The 
disc flux, the X-ray flux and the spectral distribution of illu- 
minating radiation are the same as in the lamp post model 
with rh = 0.03. 

Clearly, the approach based on local cooling cannot 
asymptotically reproduce the temperature profile given by 
the diffusion approximation. The local approach indicates 
also a certain flattening just below 10 6 K which is less im- 
portant in the titan solution. This is due to the lack of the 
true radiative transfer in this simple exercise. 

However we see that in the upper hot zone, the temper- 
ature drop occurs faster when thermal conduction is taken 




lamp post model _ 
in =0.03 



-2-10 1 
lo g( T „) 

Figure 18. The temperature profile for the lamp post model 
computed with a local approach with (solid line) and without 
(thick points) thermal conduction. All profiles are computed for 
the same parameters as the lamp post model with rh = 0.03. 
Dotted line has the same meaning as in Fig. |l]. 

into account, but finally the total optical depth of the Comp- 
ton heated zone with thermal conduction is not changed 
considerably. Therefore, we conclude that the optical depth 
of the hot layer is determined accurately enough in our ra- 
diative transfer calculations without thermal conduction. 

4.2 Comparison with other methods 

There were many computations of the vertical structure and 
of the spectrum of irradiated discs for AGN and X-ray bi- 
naries in the two past decades. The most recent ones which 
present spectra including lines from illuminated disc atmo- 
sphere with heavy elements in the case of hydrostatic equi- 
librium are NKK00 and Nayakshin & Kallman (2001), and 
Ballantyne et al. (2001). 

It is extremely difficult to compare directly our results 
to them, as our computations were performed before the 
publication of these papers, and therefore adopt other val- 
ues of the various parameters. For instance the spectral dis- 
tribution of the incident spectrum is different, a fact of most 
importance for the structure of the irradiated skin. We have 
mentioned in the introduction that the spectrum of an irra- 
diated thick layer has been computed using the three codes 
with exactly the same parameters in the case of a constant 
density medium (cf. Pequignot et al. 2001). The three com- 
puted spectra differ in their overall spectral distribution and 
in their detailed features. Concerning the model of irradi- 
ated disc in hydrostatic equilibrium, for the moment the 
only possible comparison is to stress the differences between 
the methods which are used. 

Ballantyne et al. results are limited to very high fluxes 
by their computational method, because their code has dif- 
ficulty if the hot skin has tt ^ 0.2 (see Sec. 6 Ballantyne et 
al. 2001). It precludes any comparison with their results in 
the present paper, but this will be achieved in a following 
paper, dealing with high illumination. 

NKK00 use for the thermal and ionization equilibrium 
the code XSTAR. The advantage is that the atomic physics 
is better taken into account than by us, but the drawback 
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is that the total energy balance is not correctly determined 
(cf. DAC2000 an Coupe et al. 2002). The local energy bal- 
ance is also not well determined in some cases. In particular, 
NKK00 were obliged to set the gas temperature equal to 
the effective one when it was lower. This prevents them to 
compute the temperature in the deepest irradiated layers. It 
means that above a given optical thickness, which is of the 
order of 0.2 in a case comparable to our lamp post models, 
the temperature becomes constant. This is different from 
the real temperature structure (cf. Fig. ^). Moreover, they 
do not compute like us the vertical structure of the underly- 
ing disc, but they start their integration from a height given 
by the Shakura-Sunayev (1973) vertically averaged solution. 
We show analytically in the appendix that this is valid only 
if the disc scale height is negligible with respect to the scale 
height of the irradiated skin, i.e. for low accretion rates (< 
0.01 in Eddington units). We have also a completely different 
way of handling the thermal instability, as NKK00 simply 
stop their computation when the hot solution disappears, 
while we compute the structure in the region where the cold 
solution applies. 

The same problems occur in the paper of Ballantyne et 
al. (2001), where the basis of the hot skin is fixed at the ver- 
tically averaged half-thickness of a gas-pressure dominated 
disc. Also those authors use poor description of heavy el- 
ements with only few ionization stages and excited levels. 
This treatment should underestimate the photoionization 
heating and the line cooling. 

The transfer of the continuum is performed basically in 
the same way in NKK00 and in our computations, which is 
better than the approach of Ballantyne et al. who apply the 
diffusion approximation to the transfer of the diffuse radia- 
tion (Ross & Fabian 1993). However both NKK00 and Bal- 
lantyne use the Companeets equation for the Compton diffu- 
sions, while we take the inelastic scattering into account only 
above 1 keV, through the coupling of a Monte Carlo code. 
This is clearly a disadvantage of our method. The transfer 
of the lines differ also strongly between us and them. Above 
1 keV, we use the Monte Carlo code which takes into ac- 
count Compton diffusion, and below 1 keV, we perform real 
transfer computation in the lines, but we do not take into 
account comptonization. So we have the advantage of per- 
forming real transfer in the lines, but the disadvantage of 
not taking into account Compton diffusions. This can have 
some impact on the intensities of resonance lines which are 
subject to a large number of scatterings. All other authors 
use the escape probability formalism for the lines, which 
however has the advantage to take into account the escape 
by Compton diffusions. 

In their lamp post computations, Nayakshin & Kallman 
(2001) assume a smaller X-ray flux, with the luminosity of 
the X-ray source being 20% of the integrated disc luminos- 
ity, but nevertheless it corresponds to a larger local flux 
than ours (they compute the spectrum at QRschw and not 
at WRschw)- 

NKK00 introduce the "gravity parameter" A which is 
the ratio of the vertical component of the gravity force at 
one disc scale height to the radiation pressure force (for the 
definition, see NKK00 Sect. 3.3), and they give their results 
as functions of this parameter. If we calculate the value of 
parameter A of NKK00 from our computations of the disc 
vertical structure we obtain A ~ 2 for both accretion rates. 



Indeed this parameter is not sufficient to determine the so- 
lution, one must also define the X-ray flux since A is propor- 
tional to the ratio of the viscous to the X-ray flux. So in the 
following paper, Nayakshin & Kallman (2001) present the 
NKK00 lamp post models for different accretion rates and 
X-ray fluxes. Our two models have about the same ratios of 
viscous to X-ray fluxes and about the same accretion rates 
as their models w4 and w7. Our models predict an optical 
depth of the entire hot layer of about 0.2 and 0.4, while in 
Nayakshin & Kallman the hot layer is almost twice larger. 
This can be partly due to the spectral distribution of the 
incident continuum. We adopt a power law with a slope of 
0.9 with a cut-off at 100 keV, while their power low has 
a slope of 0.8 and the continuum curves slowly above 100 
keV due to the adopted cut-off. We see clearly that a direct 
comparison is not possible. 

Finally the effect of coronal influence is studied only 
in case of purely hydrogen atmosphere in the paper of 
Kawaguchi, Shimura & Mineshige (2001). Since they do not 
take into account heavy elements, the disc does not absorb 
X-rays and produces relatively less seed photons for Comp- 
tonization. It leads to less efficient coronal cooling and pro- 
duces rather hard X-ray spectrum even if corona covers the 
whole disc. This is not the case in our calculations. The 
disc atmosphere with heavy elements has a low albedo and 
all absorbed X-ray photons are reemitted in the soft band 
as seed photons for Comptonization. Therefore, to produce 
hard X-ray spectrum the corona has to cover the disc only 
partially. 



5 CONCLUSIONS 

In this paper we have presented computations of the ver- 
tical structure and spectrum of an illuminated accretion 
disc, using a new code available for Compton thick media, 
titan/noar, and taking into account hydrostatic equilib- 
rium. The disc atmosphere is integrated vertically together 
with the disc interior. Below the illuminated atmosphere, for 
T es > 6, diffusion approximation is used. We compare two 
geometries: the lamp post model, where an X-ray source 
is situated above the disc at some distance, and the heavy 
corona model, where the corona, besides illumination, exerts 
a dynamical pressure on the disc material. We found: 

• In the case of the lamp post model, the uppermost lay- 
ers of atmosphere are heated to very high temperature which 
drops suddenly with increasing optical depth. Such a sharp 
transition is consistent with theoretical predictions and, de- 
pending on strength of illumination, may indicate thermal 
instabilities (KMT81, RC96). The radiative transfer com- 
putations prevent the temperature profile to be multivalue, 
nevertheless thermal instabilities occur as a 'wiggle' in the 
pressure profile. To avoid this problem thermal conduction 
should be taken into account and we plan to do it in the 
future work. 

• The disc scale height is important for accretion rates 
m > 0.01 (see appendix) and it influences the optical thick- 
ness of the illuminated zone. Therefore, it is essential to 
integrate carefully the cold region below the transition, in 
order to determine the fractional abundances of weakly ion- 
ized heavy elements, which influence the extreme UV and 
soft X-ray spectrum. 
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• In case of the lamp post model, the outgoing spectrum 
differs in extreme UV and soft X-ray band from the case 
of constant density slab, indicating, that it is extremely im- 
portant to estimate the temperature and the density in the 
partially ionized medium. Also K a iron line looks different 
in both cases. The absorption edge is absent contrary to the 
constant density case. 

• A heavy corona, even optically thin, completely sup- 
presses the highly ionized zone on the top of the accretion 
disc. 

• The spectrum from the disc covered by a heavy corona 
is featureless, but it depends on the expected energy distri- 
bution of incident radiation. When the X-ray spectrum is 
hard the energy balance requires that corona covers the disc 
only partially, and the spectrum from the non covered part 
becomes important. In such a case it is impossible to distin- 
guish the presence or absence of a patchy corona. When the 
X-ray spectrum is soft, we can infer the presence of a corona 
in those objects with featureless spectra. 

• Our model is in agreement with recent observations of 
K a iron line, where in some objects only the fluorescence 
component is present (NGC 5548 Yaqoob et al. 2001), in 
others, only the recombination component is present (Ark 
564 Comastri et al. 2001), or like in Mkr 205 (Reeves et al. 
2001) both components are seen. It can be explained in the 
framework of the our model by changing the illuminating 
X-ray flux. 



APPENDIX A: SEMI-ANALYTICAL RESULTS 

A crucial question is how the optical depth of X-ray heated 
zone depends on explicit calculations of the disc vertical 
structure. 

Let us consider hydrostatic equilibrium with gravity g 
and ignoring the term of radiative pressure. We can then 
simply calculate the optical depth of the hot zone as: 



hot 



Pbot^e 



(Al) 



where P(, ot is the gas pressure on the bottom of the X-ray 
heated skin, which can be expressed as: 



Pbot = - 



Fx 



(A2) 



" botC 

The ionization parameter on the bottom, E;, ot , is equal to the 
1.22(T/c/10 8 )" 3/2 (Begelman, McKee, Shields 1983), with 
the inverse Compton temperature Tic defined as usual. 

In the case of stellar atmospheres, the optical depth of 
X-ray heated layer is a function of only the hard X-ray flux 
and the shape of incident radiation (through the Tic)- 

hot FxKes 



-■'hot eg 
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In the case of an accretion disc, the gravity depends 
on the distance from the black hole and on the pressure 
scale-height of the considered layer (we assume the gravity 
to be vertically constant). Considering pressure scale-height 
of the full system i.e. the illuminated zone plus the disc 
H p = Hhot + Hd, we can simply derive: 



Fx ^e. 

E bot c Q, 2 K 



+ H d 




r/rsch— 10 

M BH =10 8 M Q - squares 
M BH =10 M G - triangles 
f = - 1 dashed line 
1=0-5 dolled line 
1 = 0.91 solid line 
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Figure Al. The ratio of scale-height of isothermal Compton 
heated skin to the disc scale-height, versus the accretion rate in 
Eddington units, for a given radius (r = lORschw) an d for differ- 
ent fractions of energy dissipated in the X-ray source: / = 0.1 - 
dashed line, / = 0.5 - dotted line, and / = 0.91 - solid line. The 
results are presented for two different masses of the black hole: 
Mb a = 10 s Mq - squares, and M BH = 1OM - triangles. Filled 
squares mark the two cases computed numerically and considered 
in detail in the previous sections m = 0.03 and m = 0.3. 



where Hhot is the scale-height of the isothermal Compton 
heated skin. If Hhot is much larger tha n H d the disc scale- 
height can be neglected. Rewriting Eq. 



A4 



with Eq. for 
constant gravity case and replacing T cor by Tic one can get: 
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For the case of variable gravity the hydrostatic equilib- 
rium should be solved together with the equation of mass 
continuity dr = K es pdz. 7v s 0< is then larger by a factor of \f2 
(the same factor as for the pressure scale-height, see Eq. |l2[) , 
and by integral of order of unity (Nayakshin 2000, Eq.2). 

If Hhot is comparable to Hd, the disc structure is im- 
portant and the optical depth of the X-ray heated skin can 
not be determined without proper disc calculations. Fig. A.1 
gives the ratio Hhot/H d versus accretion rate for the repre- 
sentative radius r = lORschw Hd is com puted by solving the 
explicit disc vertical structure (see Sec. 2.2.1) instead of as- 
suming vertically averaged Shakura & Sunayev disc (1973) 
like in NKK00. The computations are done for two black 
hole masses: 10 8 Mq, which corresponds to the AGN, and 
WMq, corresponds to GBHs. We assume different values 
of / from low illumination / = 0.1, which corresponds to 
Fx I Fdisc = 0.11, to strong illuminaton / = 0.91 which gives 
Fx/F dlsc = 10. 

One can see that the ratio is slightly higher for higher 
Fx and increases with decreasing accretion rate. Note that 
in case of GBHs the ratio is always smaller than or equal 
to unity implying that the disc structure is important for 
all accretion rates. In case of AGN for rh > 0.04 the verti- 
cal structure of the disc should be integrated down to the 
midplane to determine correctly the optical depth of the hot 
skin. 



(A4) 



Fig. A2 presents the Thomson thickness of the X-ray 
heated skin versus the accretion rate for the radius lORschw 
which is expected to be the region where the iron line is 



© 0000 RAS, MNRAS 000, 000-000 



16 A. Rozanska, A.-M. Dumont, B. Czerny, S. Collin 




001 01 0.1 1 001 01 1 I 

acc rale (In Edd. units) acc rate (in Edd. units) 



Figure A2. The dependence of optical depth of X-ray heated 
skin on the top of accretion disc on the accretion rate in Ed- 
dington units, for the given radius (r = lORschw) an d for dif- 
ferent fraction of energy dissipated in the X-ray source: / = 0.1 
- dashed line, / = 0.5 - dotted line, and / = 0.91 - solid line. 
Results are presented for two different masses of the black hole: 
Mgf{ = 10 8 Mq - squares, and Mgg = IOMq - triangles. Panel 
a represents results computed according to Nayakshin (2000) as- 
suming int egra l equal 1, and panel b is computed for constant 
gravity Eq. |A4| but including the influence of disc structure. Filled 



squares 
sidered 



mark the two solutions computed numerically and con- 
h the pteviourf sections th = 0.03 and Vh = 0.3. 



produced. We compare the analytical results of Nayakshin 
(2000), who neglect the disc scale-height and shown that this 
is in agreement with NKK00 numerical computation (see 
Nayakshin 2000 Fig.l) with ours Eq. A.4. One can see that 
taking into account proper disc calculations it is impossible 
to obtain optically thick illuminated skin even for the case 
when Fx/Fdisc = 10. Therefore, when observations imply 
that the hot layer is optically thick, it may be evidence for 
the existence of hot corona. 
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